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Abstract 

Oriented self-avoiding walks (OSAWs) on a square lattice are studied, with 
binding energies between steps that are oriented parallel across a face of the 
lattice. By means of exact enumeration and Monte Carlo simulation, we 
reconstruct the shape of the partition function and show that this system fea- 
tures a first-order phase transition from a free phase to a tight-spiral phase at 
P c = log(/i), where [i = 2.638 is the growth constant for SAWs. With Monte 
Carlo simulations we show that parallel contacts happen predominantly be- 
tween a step close to the end of the OSAW and another step nearby; this 
appears to cause the expected number of parallel contacts to saturate at large 
lengths of the OSAW. 

I. INTRODUCTION 

Many aspects of the behavior of polymers can be described by self-avoiding walks on 
a lattice. Some polymers have interactions that depend on the spatial orientation of the 
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polymer, for instance A-B polyester. Such polymers are conveniently modeled by oriented 
self-avoiding walks (OSAW) with two types of short-ranged interaction between edges de- 
pending on their relative orient at ionS^i. 

The model of investigation in this paper consists of one OSAW on a square lattice. 
Besides self- avoidance, the only interactions of the OSAW with itself occur if two steps of 
the walk are one lattice spacing apart. If the two steps have the same orientation, they are 
said to form a parallel contact, to which an energy gain of e p is attributed. If they have 
opposite orientation, they are said to form an anti-parallel contact, with an energy gain of 
e a . If (3 is the inverse temperature, and we define (3 P = —(3e p and (3 a = ~l3e a , then the 
partition function of such an oriented self-avoiding walk is given by 

Z n ((3 p ,(3 a ) = £ C n {m p ,m a )e^+^ , (1) 

m p ,m a 

where the sum is over all allowed values of the number of parallel contacts, m p , and the 
number of anti-parallel contacts, m a , and C n (m p ,m a ) is the number of configurations of 
length n with m p parallel and m a anti-parallel contacts. The limiting reduced free energy 
per step is given by 

F((3 P , (3 a ) = lim - log [Z n (p p , p a )] . (2) 

The phase diagram of this model has been studied previously!!; numerical results from 
exact series up to n = 29 edges showed the existence of three phases: a free SAW phase, 
a normal collapsed phase and a compact spiral phase. The transition from the free to the 
spiral phase was conjectured to be of first order. 

In this article we will concentrate on the case where there are only interactions between 
parallel contacts, i.e. (3 a = 0. The earlier work! rigorously proved that for this case the 
reduced limiting free energy is constant for f3 p < with value log(/x), where \x is the growth 
constant for SAW {p, = 2.638). For (3 P > the following rigorous bounds were proved: 

/3 P <F(/3 P ,0) </3 p + log(/i). (3) 
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The above results prove the existence of a phase transition for < (3 P < log(/i). Bennett- 
Wood et a% conjectured that the critical inverse temperature f3 c is near or equal the lower 
bound which, for f3 a = 0, is log(/i) ~ 1. In section [TT] we further investigate this phase 
transition by extending the exact enumeration data, by means of Monte Carlo results and 
combining them with some rigorous results on tight spirals. 

Another interesting question concerning OSAWs is the mean number of contacts. One 
of us proved that the mean number of anti-parallel contacts (m a ) ~ n in two or higher 
dimensions, where n is the number of steps of the walk!. The mean number of parallel 
contacts scales as (m p ) ~ n in three or higher dimensions, but in two dimensions the 
behavior is still an open question. Field theoretic workS predicts that in two dimensions 
(m p ) ~ log(n) in the limit n — > oo. However, Monte Carlo results for OSAWs with up 
to 3000 steps seem to indicate that {m p ) tends to a constant ~ 0.050. In section |II| we 
present the results of large-scale Monte Carlo simulations with OSAWs of up to 5000 steps, 
and investigate these results in a way that allows extrapolations to even larger n. Based on 
these results we obtain an upper bound for (m p ) in the limit n — > oo. 

II. PHASE TRANSITION TOWARDS A TIGHT SPIRAL 

Bennett-wood et ctfi enumerated all configurations up to SAWs with a length of n = 29 
and ordered them according to their number of parallel and anti-parallel contacts. We 
extended the exact enumeration of the OSAWs with parallel contacts, and obtained all 
values for C n (m p ), the number of OSAWs consisting of n steps and having m p parallel 
contacts, up to n = 34. 

In our enumeration program, we start with generating all OSAWs of length / < n with 
a parallel contact between the first and the last step. For each walk w, we determine the 
number of parallel contacts m w . We also determine Mi(w,ti,rrii), the number of extensions 
on the inside end of walk w with length £j < n — I and m,; parallel contacts with either 
itself or w, and M (w,t ,m ), the corresponding quantity for the extensions on the outside 
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end. Since the walk w prevents contacts between the inside and outside extensions, the total 
number of OSAWs of length n with m p parallel contacts is given by 

Z n (m p ) = —^2 H Mi(w,ti,mi) M (w,t ,m ) 5(m p ,m w + m, + m ). (4) 

m P w l +ti+ t =n 

The prefactor in this equation corrects for the fact that there are m p different walks w from 
which we can generate the same OSAW with m p parallel contacts. Exploiting rotational 
and mirror symmetry, we enumerated all OSAWs of length n < 34 with one or more parallel 
contacts in a run of about two weeks on a four-processor DEC alpha workstation. Finally, 
the number of OSAWs without parallel contacts is obtained by subtraction from the total 
number of OSAWs (from Ref. 7). 

The results are given in table |, and plotted as the solid lines in figure [TJ, where 
log(C n (m p )) is plotted as a function of m p . The figure shows that up to n = 34, the 
number of configurations as a function of the number of parallel contacts first drops quickly 
with a factor p n , but then, over the whole range 1 < m p < m max , falls off exponentially with 
the same exponent q n . The partition function Z n (f3 p ) is thus described well by 

C n (l) = Pn C n (0) 

C n {m) « C n {\) ■ exp(-g n (m - 1)), (5) 

where p n and q n are n— dependent parameters. 

To extend the graph presenting the partition function beyond n = 34 by means of exact 
enumeration is very hard. However, the left part of this graph for much larger n can be 
obtained statistically by means of Monte Carlo simulations: OSAWs are randomly generated 
with the pivot algorithm!, and a histogram is made of the number of parallel contacts of 
these OSAWs. This gives us a direct measurement of C n (m p )/Z n (0) for a small number of 
parallel contacts. In our Monte Carlo simulations, we thermalized over 10 7 pivot moves, 
followed by 10 8 moves to gather statistics; statistical errors were obtained by repeating the 
whole procedure 10 times. The results are shown in table [TT]; the density of OSAWs with 
more than ~ 10 parallel contacts is so small that they will most likely never be generated, 
and we only obtain an upper bound for them. A good approximation for Z n (0) is known: 

4 



Z n (0) « {A/A^V 1 (6) 

where fi = 2.638, 7 S = 43/321, and A = 1.7710. The factor of a fourth is due to the fact 
that we count OSAWs that are equivalent after rotation only once. The Monte Carlo results 
from table [TJ for n = 50, 60, 70, 80, 90 and 100, multiplied by Z n (0), are plotted as circles in 
the left side of figure [TJ. 

Also the utmost right point of the graph can be obtained, as there the only relevant 
configurations are tight spirals. The corners of a tight spiral are reached after n = k, k + 
1, 2k + 2, 2k + 4, 3k + 6, 3k + 9, ■ • ■ steps, i.e., at n = ik + i 2 or n = ik + i(i — 1), where k 
is the number of steps in the same direction at the inner end of the tight spiral, and % is a 
positive integer. Each additional step of the tight spiral adds one parallel contact, except 
steps before and after a corner. Thus, the number of parallel contacts m max for an OSAW 
of length n is given by 



if (n < 2k) : m r 
if (n > 2k) : m r , 



n-2k + 3 



k 2 k 

! n H 

4 2 



k 2 , k 

'n H 1 

4 2 



'n + 



{k-1) 2 (k-1) 



'n + 



(k-V 



(k-1) 



(7) 



4 2 

where square brackets denote the Entier function ([x] is the largest integer not larger than x). 
The number of parallel contacts of a 'rectangular' tight spiral (with k > 1) does never exceed 
that of the 'square' tight spiral (with k = 1), but can be equal, adding to the degeneracy 
of the ground state. Additional ground states can be generated by removing steps from 
the inside and adding them to the outside end, until the corner is reached. Also, if the 
tight spiral ends at a corner or one or two steps further, additional groundstates arise by 
rearranging these last steps. 

We enumerated all OSAWs with m max parallel contacts and length up to n = 50, and 
confirmed that all groundstates can be generated with these operations. Assuming that no 
new types of degenerate groundstates arise after n = 50, we calculated the degeneracy of the 
ground state for lengths up to a million steps, and observed that the degeneracy fluctuates 



between 4 (for a complete 'square' tight spiral) and c m n 3 / 4 with c m = 5.3, whereas the 
expected degeneracy grows as c a n 3 ^ with c a = 2.1. For n = 50, 60, 70, 80, 90 and 100, there 
are 140, 40, 16, 4, 16, and 8 configurations with the maximum number of parallel contacts. 
We have added these results of the tight spirals in figure [l] as squares. 

For n 3> 34, the Monte Carlo data in figure [1] for small m p do not extrapolate to the 
exact results for tight spirals, but point below, which suggest that eq. (|) is an upper bound 
for n 3> 34. The dotted lines in figure [1] represent these upper bounds. We cannot exclude 
the possibility that for n ^> 34 the partition function initially stays below these dotted lines, 
then increases and crosses this dotted line for intermediate values of m p , and finally reaches 
the exact result for tight spirals; however, we think that that scenario is unlikely, and the 
results concerning long OSAWs in the remainder of this section are based on the assumption 
that the dotted lines in figure p] represent upper bounds. 

For n < 34 we know C n (0) and C n (l) by exact enumeration, and for n = 50, 60, 70, 80, 
90, 100, 1000 and 2000 we know C n {0)/Z n and C n {l)/Z n accurately from the Monte Carlo 
simulations. This enables us to compute p n in eq. (|5|) for all these values of n. For large n, 
p n converges to a constant value around 0.031. To extract the specific heat and density of 
parallel contacts, we used a fit to p n which is given by 

Pn-Poo^ (1/Vn), (8) 

where p^ = 0.031 ± 0.002. We can obtain the values q n in eq. |^ from equations (U), (fTf), 
and (H), as 

_ log(C n (l)) - \og{C n {m max )) ^ \og{Z n ) + log(p n ) - 3/4 log(w) 

For n=1000 and 2000, this equation predicts that g n =1.099 and 1.060, respectively, whereas 
the Monte Carlo results in table [IT] for C n (l)/C n (5) indicate that the slope of log(C n (m p )) 
corresponds to values of q n ~ 1.4; for larger values of n the curves of log(C n (m p )) versus 
m p initially point below the point corresponding to the tight-spiral configuration, and thus 
must bend upwards at larger m p . 
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For n up to 34 we plotted in figure § the specific heat, defined by ftx = —d 2 F/d(3 2 , and 
in figure |3| we plotted the density of parallel contacts (m p )/n, as a function of the inverse 
temperature [3. In both figures we added the graphs for n = 50, 100, 200, 500, 1000, 2000 
and 10,000, obtained from eq. ([5]), as dotted lines. In figure Q, the value of (3 where the 
peak of the specific heat is located is moving backward to (3 = log(/x), as is the point where 
{m p )/n is increasing steeply in figure ||. The jump in the density of parallel contacts (i.e., 
the energy density) is increasing with increasing n, indicating a first order transition. In 
fact, assuming eq. @ one can show analytically that in the limit n — > oo the function 
(m p )/n approaches the Heaviside stepfunction 0(log(/i)), and this still holds if eq. (||) is an 
upper bound rather than an exact expression in the regime between tight spirals and walks 
with few parallel contacts. Both the specific heat and the density of parallel contacts are 
insensitive to the fact mentioned earlier, that the curve starts somewhat steeper at small 
m p and thus must bend up at larger m p . If anything, they will increase the peak value of 
the specific heat, and the steepness of the density curve. 

Another way to estimate the transition point is to look at the zeroes of the partition 
function!'!. The partition function of an OSAW of n steps with m p parallel contacts is 
a polynomial of degree m max (the maximum number of parallel contacts) in the variable 
x = e@, hence it can be conveniently written in terms of its n roots r mp in the complex 
plane: 

Z n (x) = C n (0) J] (1 - (x/r mp )) (10) 

m p =l 

and the free energy per steps 

j 1 ITlmax 

F n (x) = - log(C„(0)) + Ml - (*Av»J). (11) 

The coefficients C n (m) are real and non-negative, hence none of the roots lies on the real 
positive axis, but for n —>■ oo they will cross it at some point x c < /i, since we rigorously 
know the existence of a phase transition. 

We calculated the zeroes of the partition function corresponding to the exact data up 
to n = 34 and they are plotted in figure ^a. The roots seem to lie in nearly perfect circles 
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for every n, but the radius decreases with increasing n. The n roots nearest to the real 
positive axis approach the real axis along a nearly straight line. In figure f|b, we plotted the 
real part of the root nearest to the real axis for n = 25. .34, against 1/n. Again, the figures 
are consistent with a transition at x c ~ 2.5. 

III. NUMBER OF PARALLEL CONTACTS FOR (3 = 

The second major topic of this paper is to investigate the behavior of the number of 
parallel contacts m p in the limit n — > oo. In figure [5] we have plotted the behavior of (m p ) 
as a function of n, obtained from eq. (|5]), which we proposed to be an upper bound. The 
upper bound reaches asymptotically the value m p = 0.08. Clearly, the earlier mentioned 
fact that the curve has a somewhat steeper slope at large n and small m p has impact on 
(m p ), as these configurations are dominant at /3 = 0. Therefore we do not use eq. (|5|) in the 
remainder of this section. With Monte Carlo simulations we have determined the expected 
number of parallel contacts (m p ) as a function of n. The results are given in table [TTT] and 
figure |5|, and are in agreement with results published earlier by one of usi, but extend to 
larger values of n. The Monte Carlo results seem to converge to a value around 0.05. 

To understand the underlying physics in the regime j3 = better, we took a closer look 
on where the parallel contacts are made, and relate this to other types of SAWs. Consider 
an oriented OSAW of length n, with a parallel contact between the steps i and j of the walk. 
The sequence of steps from i to j constitute a polygon of length I — j — i + 1, if one of the 
two steps that form a contact is rotated 90 degrees to close the polygon. The remaining 
sequences of steps from to i and from j to n are two self-avoiding walks of length % and 
n — j, respectively. These two SAWs can be combined into one self-avoiding two-legged star: 
a SAW of length n — I, on which one special point (the origin of the two-legged star) is 
marked. Note that, since the two SAWs are separated by the loop, one being located on the 
inside of the loop and one on the outside, the two-legged star is always self-avoiding. The 
mapping of an OSAW with one parallel contact into a rooted polygon plus a two-legged star 
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is illustrated in figure ||. 

If an OSAW has more than one parallel contacts, then we can map this OSAW onto 
different combinations of a rooted polygon plus a two-legged star. In general, if the OSAW 
has m p parallel contacts, there are m p such mappings into a rooted polygon plus a two- 
legged star. The reverse mapping, i.e. the combination of a two-legged star plus a rooted 
polygon into an OSAW with a parallel contact, is not guaranteed to result in an OSAW with 
a parallel contact, as they might cross. Therefore, the total number of rooted polygons of 
length / times the total number of two-legged stars of length n — I, summed over all I, is an 
upper bound to the number of OSAWs of length n, multiplied by the expectation value of 
the number of parallel contacts for these walks. 

Let us define f(n, I) as the probability that a two-legged star of length n — I if combined 
with a rooted polygon of length I, results in an OSAW. Then we can write 

(m p )Z n = m p C n (m p ) = ^ PiS n ^if(n, I) (12) 

mp l 

where Z n , P n and S n are the number of OSAWs, rooted polygons and two-legged stars of 
length n, respectively. 



We know that, for large n: 



Z n « fxV- 1 (13) 
S n « // n n 7s (14) 
P n « fi n n a - 2 (15) 



Combining this with fll2|) leads to: 

l a - 2 {n-l)^f{n,l) 



K) = E ,,'V ' (16) 

We can obtain insight in the behavior of the function f(n,l) by means of Monte Carlo 
simulations. OSAWs are sampled randomly, and for each parallel contact the loop length 
I = \j ~ i + 1| is determined, where % and j are the steps making the parallel contact. This 
procedure gives us (m p )(l), the expectation value of the number of parallel contacts with 
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loop length /. Results for OSAWs with a length of n =200, 500, 1000, 2000, and 5000 are 
plotted in figure [7]. (m p )(l) shows a power-law behavior, where the length n of the OSAW is 
an upper bound to the length I of the loop. Important however is that, besides this obvious 
dependence, the total length n does not appear to have any influence on the behavior of 
(m p )(l), and this quantity is well described by a power-law: 



(m p )(l) « k r a ' (17) 



Numerically, we find: 



k = 0.35 ± 0.1 (18) 
ai = 1.65 ± 0.05 (19) 

To obtain the mean number of parallel contacts (m p ) we sum over all possible (even) lengths 
I of the rooted polygon: 

n 

K> = EK>(0«* E rQi - (2°) 

I 1=8 

For n — > oo the right hand side equals a constant times the function ((ai), which converges 
to a constant for a\ > 1. This implies again that (m p ) tends to a constant in agreement 
with earlier Monte Carlo results of Flesiai. 

The fact that (m p ) is constant implies that the SAW critical exponent 7 is constant in 
the free and repulsive regime (i.e. for (3 < 0), and presumably until the transition. For the 
exponent 7 to change with (3, the exponent a\ should be < 1, since this will cause the ( 
function to diverge, but this is not supported by our numerical results in Fig. 7. 

It is possible to put an upper limit to how far {m p ) will still increase if n is increased 
above 5000: figure [7] shows that the contribution of loops with a length below 1000 certainly 
has converged for n = 5000, thus (m p )(oo) - (m p )(5000) < k ■ EL1000 ^ < 10 ~ 5 - 

A different approach which estimates both the number of parallel and anti-parallel con- 
tacts is to use the similarity between an OSAW and a twin-tailed tadpole. Consider an 
OSAW with a contact between steps i and j of the walk. If we add a new edge between 
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steps % and j we obtain an object which we will call a twin-tailed loop (see figure 8). A 
twin-tailed loop differs from a non-uniform twin-tailed tadpole only by one edge, and has 
the same asymptotic behavior. If the contact is parallel then the twin tailed loop has one 
tail inside the loop and the other outside (see fig 8a), while if the contact is anti-parallel 
both tails are outside (see fig 8b). This is of course only true in two dimensions. Each 
OSAW with m contacts can be mapped into m distinct twin-tailed loops. If T n is the total 
number of twin-tailed loops of total length n then it follows that 

^ = ^^0^171)). (21) 

Dividing both sides by Z n , where Z n is the partition function of SAWs, it follows that 

(m) = T n /Z n (22) 

Asymptotically, Z n ~ /i n ?2 7s_1 , where 7 S is the exponent for SAWs. Lookmann0 proved 
that twin-tailed tadpoles have the same growth constant fi as SAWs and that the exponent 
7 is 7 = 7 S + 1. The same kind of proof holds for twin-tailed loops. Replacing these results 
in eq. ( |2"2"D implies the known result (m) ~ n. 

Consider now the parallel and the anti-parallel case separately. Twin-tailed loops with 
both tails outside the loop are the dominant configurations, so they have the exponent 7 of 
the total set, i.e. 7 = 7 S + 1. This implies as previously that (m a ) ~ n as was proven by 
one of us!. 

Parallel contacts correspond to the subset T* of twin-tailed loops with one tail on the 
inside and one on the outside of the loop. The question is, what is the value of the exponent 
7 (let us call this exponent j t ) for this subset T*l Simple tadpoles {i.e. tadpoles with only 
one tail) have the same 7 as SAWs0. Since one element of T* can be constructed from a 
simple tadpole by adding one edge inside the loop, it follows that j t > j s . On the other 
hand, since T* is a subset of the set of twin-tailed loops, it follows that 7t < 7 S + 1, and this 
inequality can be made strict by considering that (m p ) ~ o(n), see Bennett-Wood et all 

We can gain insight in this matter by randomly generating OSAWs of length n, and 
for each parallel contact determining the length t of the inside tail. Note that if a parallel 
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contact is formed between steps % and j of the OSAW, the steps from % to j form a loop, 
and 'inside' and 'outside' tails refer to inside or outside this loop. The results are plotted in 
figure 9. Extrapolating these results we estimate that the fraction of twin-tailed loops with 
length t of the inside tail is decreasing as 

(rn p )(*)~fc t r at (23) 

where at = 1.6 ± 0.1. The parameters a\ and at are within each others statistical errors 
and are probably the same. As in Eq. ( p3|) the parameter at exceeds 1, J2t{ m p{t)) wm n °t 
be more than a constant times m p (t = 0). This implies that T* asymptotically seems to 
behave as simple tadpoles which have the same 7 as SAWs. If we assume, based on these 
numerical results and intuitive arguments, that the twin-tailed loops with one tail inside and 
one outside behave as simple tadpoles then 7t = 7 S , which would imply that (m p ) approaches 
a constant. 
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TABLES 
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TABLE II. Monte Car] 


lo results for the i 


density of OSAWs 


of length n with 


m p parallel contacts. 


n 


m p =0 


1 


2 


3 


4 


50 


0.97763(1) 


0.01841(1) 


0.003209(8) 


0.000599(3) 


0.000120(1) 


60 


0.97603(2) 


0.01954(2) 


0.003555(7) 


0.000696(3) 


0.0001426(8) 


70 


0.97479(3) 


0.02039(3) 


0.003832(7) 


0.000780(4) 


0.000164(1) 


80 


0.97368(2) 


0.02118(2) 


0.004067(10) 


0.000840(6) 
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0.0000052(3) 
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0.0000300(5) 


0.0000062(2) 


0.0000013(1) 


0.00000026(5) 
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70 


0.0000346(6) 


0.0000080(3) 


0.0000015(2) 


0.00000021(4) 


0.00000009(4) 


80 


0.0000401(8) 


0.0000087(5) 


0.0000019(2) 


0.00000032(7) 


0.00000014(3) 


90 


0.000044(1) 


0.0000109(3) 


0.0000020(2) 


0.00000045(6) 


0.00000011(4) 


100 


0.000047(1) 


0.0000112(5) 


0.0000029(2) 


0.0000007(1) 


0.00000012(4) 


1000 


0.000102(8) 


0.000032(6) 


0.000009(2) 


0.0000007(4) 




2000 


0.00012(2) 


0.000024(5) 


0.000009(3) 


0.0000010(4) 
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TABLE III. 


Monte Carlo data for (m p ) 


the expected number of to 


;al parallel contacts. 


n 


(m p ) 


n 


(m p ) 


9 


0.001966(3) 


10 


0.00505(3) 


11 


0.00450(5) 


12 


0.00715(3) 


13 


0.00698(2) 


14 


0.00918(3) 


15 


0.00921(3) 


16 


0.01106(4) 


17 


0.01118(3) 


18 


0.01274(4) 


19 


0.01293(2) 


20 


0.01429(3) 


21 


0.01446(4) 


22 


0.01577(4) 


23 


0.01592(2) 


24 


0.01693(7) 


28 


0.01925(2) 


29 


0.01953(7) 


30 


0.02025(3) 


38 


0.02358(5) 


39 


0.02389(7) 


40 


0.02431(8) 


41 


0.02448(4) 


48 


0.02667(5) 


49 


0.02690(8) 


50 


0.02731(7) 


70 


0.03129(7) 


71 


0.03141(5) 


80 


0.03281(7) 


90 


0.0338(4) 


99 


0.0350(5) 


120 


0.0372(4) 


150 


0.0385(4) 


200 


0.0406(4) 


300 


0.0429(4) 


400 


0.0446(4) 


500 


0.0462(7) 


700 


0.0471(6) 


1000 


0.0492(9) 


1500 


0.0493(8) 


2000 


0.0497(8) 


3000 


0.0497(9) 


5000 


0.0514(3) 
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FIGURES 

FIG. 1. A graphical representation of the partition function: the logarithm of the number 
of OSAWs is plotted as a function of its number of parallel contacts. Solid lines are data for 
n = 11. .34, obtained from exact enumeration, circles are data for n = 50, 60, 70, 80, 90, and 100, 
obtained from Monte Carlo simulations, squares from properties of tight spirals, and the dotted 
lines are connecting the Monte Carlo results with the corresponding results for tight spirals. 

FIG. 2. Specific heat as a function of inverse temperature (3. In the direction of increasing peak 
value, the curves are obtained for n=25, 30, and 34 from exact enumeration (solid lines) and for 
n=50, 100, 200, 500, 1000, 2000 and 10000 from eq. (4) (dashed lines). 

FIG. 3. density of parallel contacts, as a function of inverse temperature (3. In the direction 
of increasing density, the curves are obtained for n=25, 30, and 34 from exact enumeration (solid 
lines) and for n=50, 100, 200, 500, 1000, 2000 and 10000 from eq. (4) (dashed lines). 

FIG. 4. Left: zeroes of the polynomial of the partition function for n = 25.. 34. Right: The 
zeroes of the n-th root approach the real axis nearly along a straight line, crossing the real axis at 

X c R2 2.5. 

FIG. 5. Expected number of parallel contacts, as a function of length. The circles with error 
bars are Monte Carlo measurements, the solid line results from eq. (4) and is an upper bound. 

FIG. 6. Decomposition of an OSAW into a loop and a two-legged star 

FIG. 7. probability that an OSAW has a parallel contact with a loop of length I, for OSAWs 
with a total length n=500, 1000, 2000 and 5000. For each parallel contact, the loop length I is 
defined as I = \ j — i + 1|, where i and j are the steps of the OSAW making a parallel contact. 
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FIG. 8. An OSAW with a contact can be transformed into a twin-tailed loop by adding one 
step. If the contact was parallel, the twin-tailed loop has one tail on the inside and one on the 
outside of the loop (see figure a). If the contact was anti-parallel, both tails are located on the 
outside of the loop (see figure b) . 

FIG. 9. probability that an OSAW has a parallel contact with an inside tail of length t, for 
OSAWs with a total length w=500, 1000, 2000 and 5000. For each parallel contact, the steps i up 
to j form a loop, where i and j are the steps of the OSAW making a parallel contact. The inside 
tail is defined as those steps of the OSAW that are located within this loop. 
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